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We discuss the recently observed 'excesses' in cosmic ray electron and positron 
fluxes which have been widely interpreted as signals of dark matter. By consid- 
ering the production and acceleration of secondary electrons and positrons in 
nearby supernova remnants, we predict an additional, harder component that 
becomes dominant at high energies. The unknown spatial distribution of the 
supernova remnants introduces a stochastic uncertainty which we estimate an- 
alytically. Fitting the prediction for different source distributions to the total 
electron + positron flux measured by Fermi-LAT fixes all free parameters and 
allows us to 'postdict' the rise in the positron fraction seen by PAMELA. A 
similar rise in the B/C ratio is predicted at high energies. 
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1. Introduction 

The recent findings of a rise in the positron fraction by PAMELA 1 and an 
excess in the total e + + e~ flux by Fermi-LAT 2 have generated a lot of in- 
terest because these might be signals of dark matter annihilation or decay. 
However the expected signals are much smaller than those observed and 
models which attempt to circumvent this by invoking exotic new physics 
are increasingly constrained by 7-ray and radio observations. It is thus im- 
portant to investigate more prosaic astrophysical explanations, e.g. nearby 
pulsars 3 ' 4 or supernova remnants (SNRs). In this context we wish to revisit 
the predictions of the standard SNR origin model of galactic cosmic rays 
(GCRs) wherein only primary electrons and nuclei are accelerated in SNRs 
by diffusive shock acceleration (DSA). 5 

We consider the recent proposal 6 that some of the secondary e + made 
in the SNRs are also accelerated by DSA. By accounting for the discrete 
distribution of the SNRs, we are then able to match the predicted fluxes to 
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the measured e + + e~ flux. Having thus fixed all the free parameters of our 
model we can make an independent 'postdiction' for the positron fraction 
which agrees remarkably well with the PAMELA data. We also predict a 
similar rise for the B/C ratio as a definitive test of our model. 

2. The discreteness of sources 

Although SNRs can grow to become quite large (over <~ 100 pc in diameter) 
towards the end of their lifetime, they are still much smaller than the Kpc 
scales relevant for GCR propagation. Furthermore, the bulk of the particles 
accelerated by a SNR will be released towards the end of the Sedov-Taylor 
phase when they can no longer be confined by magnetic turbulence. We are 
therefore not dealing with a continuous distribution of sources within the 
galactic disk but with a large number of discrete sources, both in space and 
time. However, as we do not know the exact distribution, this discreteness 
introduces an uncertainty in the predicted total e + + e~ flux. However, a 
generalised central limit theorem and the analytical form of the Green's 
function for the diffusion-energy loss problem can be used to infer the 
statistical properties of the flux, i.e. its average and quantilcs. Our results 
arc in accordance with earlier studies using the Monte Carlo approach. 7 ~ 9 
The flux from a source i that injected e~ or e + a time t{ ago at a distance 
Li from the observer is given by the Green's function 10 G<n s k(E, Li, t{) of 
the diffusion equation: J t {E) = c/ (4ir)G disk (E, L u ti). The flux J(E) of A 
identical sources at distances {Li} and times {ti} is just the sum of the 
individual fluxes, 



As a function of the random variables L and t, the Green's function 
Z = Gdisk(E, L, t) is itself a random variable with probability density fz, 
expectation value \iz and standard deviation <Jz- If the central limit the- 
orem were applicable, the fluxes J for different realisations of the same 
source density at a fixed energy E would follow a normal distribution with 
mean [ij = c/(4tt)N ' \iz and standard deviation aj — c/(4n)^/Naz ■ 

In the case under consideration, L and t are assumed to be independent 
random variables with probability densities fi, and ft, respectively. We 
choose the sources to be distributed homogeneously in a ring with inner and 
outer radius n and r 2 , so Jl = 2L/ (rf— r^) for r\ < r < r 2 and otherwise. 
For simplicity we also assume the source rate to be constant in time, i.e. 
ft = 1 /imax for < t < i max and otherwise, up to a maximum time set by 



N 



N 




(1) 
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the minimum energy E m - m under consideration: i max = (fro-^min) -1 - Here, 
6 is the normalisation of the energy loss rate of the e~/e + : b(E) = boE 2 . 

For E > 10 GeV and n — > (the distance to the nearest source is not 
limited physically) we find the expectation value: 

n T - C , 1 A_g2.£r7-l+(l-g)/2 F (^) (2) 

4tt 6o(l - S) *r 2 2 t max r ^-,5-1 ^ 

where the interstellar diffusion coefficient Do = 10 28 cm 2 s -1 at 73 = 1 GeV, 
and (5 ~ 0.6 is its spectral index (D cx E 5 ). We set b = 10~ 16 GeV" 1 s" 1 
and T2 = 15kpc as is standard. The total number TV of sources that are 
needed to reproduce the (observed) number N ~ 300 of SNRs active in 
the Galaxy at any given time depends on the average lifetime of a SNR, 
t snr, which is suggested to be <~ 10 4 yr, 11 hence TV = A/"t max /r SNR = 10 6 . 
With Q = 8.4 x 10 49 GeV -1 , the expectation value, 

E 3 hj ~ 150GeV _1 cm- 2 s _1 sr _1 , (3) 

closely matches the featureless E~ 3 spectrum measured by Fermi-LAT. 2 

Whereas the expectation value [ij is well-defined, the variance ct 2 di- 
verges because fz{z) has a long power law tail. Nevertheless, instead of the 
central limit theorem, one can apply a generalised version 12 which allows 
determination of the distribution of J from the asymptotic behaviour of 
fz{z) for z — > oo. The distribution of J is a so-called stable distribution 13 
— an asymmetric generalisation of a Gaussian which also exhibits a power 
law behaviour for large values of J. The uncertainty interval around the 
expectation value [ij = c/(4Tr)Nfiz can then be defined by quantiles of the 
stable distribution. The energy dependence of the quantiles turns out to 
be always harder than that of the expectation value, so the fluctuations 
are growing with energy. Figure 1 shows the uncertainty bands around the 
expectation value for the flux fij, together with the expected fluxes from 
50 random realisations of the above source distribution. 

To calculate the fluxes of e~ and e + at Earth, we perform a Monte Carlo 
calculation by considering a large number of realisations of randomly dis- 
tributed sources according to a probability density function which reflects 
our present knowledge concerning the distribution of SNRs in the Galaxy. 
We emphasise that the better the expected flux of e~ and e+ from such a 
realisation of the source density matches the measured fluxes, the closer is 
the underlying distribution of sources likely to be to the actual one. We do 
not consider any variations between the SNRs but assume a prototypical 
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Fig. 1. Fluxes of high energy e~ from ensembles of sources uniformly distributed in a 
disk around the observer. The solid line denotes the expectation value for the sum of 
fluxes from N discrete, transient sources while the dashed lines show the expectation 
values if the sources are limited to a ring with inner radius r\ (normalised to the ex- 
pectation value for n = at lGeV). The coloured bands quantify the fluctuations and 
contain, respectively, 68 % and 90 % of the calculated fluxes. The fluxes from 50 random 
realisations of an ensemble of N individual sources are shown by the thin grey lines. 

set of source parameters which we determine from a compilation of 7-ray 
SNRs, see Sec. 3. Of course all SNRs are not the same, however variations 
of the source parameters would only introduce additional fluctuations in 
the fluxes without altering their average. 

For a more realistic distribution of source distances, we model the SNR 
density by a logarithmic spiral tracing the gas density (see Fig. 2) and 
also include the known radial distribution of SNRs in the Galaxy from 
radio surveys. We transform to a coordinate system centred on the Sun and 
integrate over azimuth such that our calculated probability density (Fig. 3) 
encodes the average surface density of SNRs at a particular distance. 

3. Primary electrons 

The injection of primary e~ by SNRs is parametrised as: 

*.-=«--(cJv)~ Vw " < 4) 

The spectral index 7 and the cut-off energy E cut can be obtained from the 
spectral indices of SNRs in 7-rays as measured by Imaging Air Cerenkov 
Telescopes (IACTs), like HESS, MAGIC and VERITAS. (We assume the 
same spectral index for the hadronic and the electronic components, as pre- 
dicted for DSA.) We have compiled 14 a list of all SNRs detected by IACTs 
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Fig. 2. The adopted distribution of Fig. 3. The probability density for the 
SNRs in the Galaxy — the cross denotes distance of a SNR from the Sun. 
the position of the Sun. 
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Fig. 4. Predicted spectra of e~ and e + with data from Fermi— LAT 2 (red circles) and 
HESS 17,18 (blue squares & green triangles). Left: Primary e~ after propagation to Earth. 
Right: Secondary e~ and e + from cosmic ray interactions, created during propagation 
(dashed line) and created during acceleration in SNRs (full lines). 

and find the typical values to be 7 = 2.4 and E cut ~ 20TeV. The normal- 
isation R°_ is determined by fitting the electron flux at Earth from our 
Monte Carlo computation to the preliminary measurement by PAMELA at 
10 GeV; 15 the secondary fluxes can be neglected for this purpose. We find 
R° = 1.8 x 10 50 GeV^ 1 for 7 = 2.4 which corresponds to a total injection 
energy of / dE E R e - (E) ~ 7 x 10 47 erg. Solar modulation is treated in the 
spherical force-field approximation 16 with a potential of <f> = 600 MV. 

The primary electron fluxes on Earth from a large number of realisa- 
tions of the source distribution are shown in the left panel of Fig. 4 and 
clearly exhibit a deficit with respect to the measurements by Fermi-LAT 
and HESS. 
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4. Additional positrons from the acceleration of secondaries 

It is usually assumed that secondary e~ and e + are produced only through 
spallation (mainly on interstellar hydrogen and helium) during the galactic 
propagation of hadronic cosmic rays. At energies above 10 GcV, the spec- 
trum of such secondaries after propagation should be softer than that of 
the primaries (by the power 5, the spectral index of the diffusion coeffi- 
cient), and the positron fraction should therefore be falling. We calculate 
secondary e + from propagation of protons and nuclei following Rcf 19, 
using however the propagation parameters specified above. 

It was recently suggested 6 that the acceleration of secondary e~ and 
e+ in the cosmic ray sources, i.e. SNRs, could lead to a harder source 
spectrum of e + that can explain the observed rise in the positron fraction. 
As the volume in which particles that participate in DSA grows with the 
diffusion coefficient -Deohm (p) <x P, the fraction of the secondaries which can 
participate in DSA also increases linearly with momentum. The resulting 
source spectrum R± of secondary e~ and e + is then a sum of two power 
laws, corresponding to the unaccelerated and the accelerated secondaries: 



R± ~ R°± p~ 



7 


1 + 


( ' )1 






\ Across / . 



(5) 



D(p CIOSS ) = - — ; 0U1 )t ; . (6) 



where the "cross-over" momentum, p C ross, satisfies 

3 rMfT SNR 
4( 7 + 2)(l/c! + r 2 

As has been noted, 6 this mechanism is most efficient for old SNRs where 
field amplification by the shock wave is not very effective anymore. We 
therfore introduce a fudge factor Kb that parameterises the effect of the 
smaller field amplification on the otherwise Bohm-like diffusion coefficient 
in the SNR, 

D {E) =3,xlft, (-§) "'(A).-,-. (7, 

A break in the source spectrum occurs at pt>reak because the growth of 
the acceleration zone is bounded by the physics size of the SNR. The source 
spectrum R± thus returns to a p~ 7 dependence around p = pbroak- At even 
higher energies the secondary spectrum cuts off at the same E cut as the 
primary e~ (see Sec. 3). 

Following Refs. 6,20, the parameters are chosen to be: U\ — 0.5 x 
10 s cms -1 , n gas ,i = 2 cm~ 3 , B = 1 [iG. Instead of fixing the normalisation 
of the injection spectrum ad hoc, we determine it from 7-ray observations 
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Fig. 5. Left: Predicted spectra of the total flux of primary and secondary e~ and e + . 
Right: The predicted positron fraction after fixing all free parameters using the total 
(e + + e~) flux, compared with data from PAMELA. 1 



of SNRs. Knowing the cross-sections for the production of 7-rays (from the 
decay of neutral pions) and the production of secondary e~ and e + (from 
the decay of charged pions), we determine the normalisation 

R°, = 7.4 x 10 48 ( ( ) Gey- 1 . (8) 

+ ViOV/ \5.7x lO^s-iTeV- 1 / KJ 

where Q° is the observed typical luminosity of SNR in 7-rays as determined 
from the same compilation of data as above. 

The right panel of Fig. 4 shows the expected flux of secondary e~ and 
e + for 30 realizations of the possible distribution of SNRs in our Galaxy. 
Clearly this component can potentially match the high energy Fermi-LAT 
and HESS data. We also show the secondary e~ and e + from the propaga- 
tion of protons and nuclei, which is subdominant. 



5. Total electron positron flux and positron fraction 

Figure 5 shows the total (e + + e~) flux obtained by adding the primary 
e~ and the secondary e~ and e + (see. Fig. 4). As explained above, we 
selected just those (3) realisations of the source distribution which give the 
best fit to the measurements by Fermi-LAT and HESS. By fitting to the 
(e + + e - ) flux we have also fixed the value for the diffusion coefficient near 
the SNR shock wave (which determines the ratio of the accelerated to the 
unaccelerated secondaries) to about 15 times the Bohm value i.e. Kb — 15. 

Having thus fixed all free parameters we can make an independent 
prediction for the positron fraction which is shown in the right panel of 
Fig. 5. It is in very good agreement with the PAMELA measurements 
above > 10 GcV; we do not expect to match the data below 10 GeV be- 
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Fig. 6. The predicted B/C ratio in cosmic rays (from Ref. 22) for the 'leaky box model' 
with production of secondaries during propagation only (dot-dashed line), and including 
production and acceleration of secondaries in nearby SNRs (solid lines) for values of the 
diffusion coefficient near the shock wave which best fit the spectrum (see Sec. 5). The 
dashed line corresponds to the value of the diffusion co-efficient required to fit ATIC-2 
data on Ti/Fe (from Ref. 22), along with its Icr and 2<r error bands. The data points are 
from HEAO-3-C2 (circles), 24 ATIC-2 (triangles) 25 and CREAM (squares). 26 

cause our analytic approach neglects complications such as convection and 
reacceleration in the ISM which can be important at these energies. 

6. Nuclear secondary-to-primary ratios 

The acceleration of secondaries described above modifes not only the spec- 
tra of secondary e~ & e + but also of other charged secondaries. This pro- 
vides a direct test of this explanation for the rise in the positron fraction, e.g. 
a rise in the antiproton-to-proton ratio is also predicted 20 (consistent with 
the measurements by PAMELA 21 so far). We have argued 22 that nuclear 
secondary-to-primary ratios can be used to not only test this model but also 
to discriminate against alternative explanations such as dark matter anni- 
hilation/decay or nearby pulsars which cannot affect nuclear abundances. 

We have calculated the rise in nuclear secondary-to-primary ratios ex- 
pected from the acceleration of secondaries, in particular, we have used 
recent ATIC-2 data on the titanium-to-iron ratio which exhibit a rise 
around 100 GeV/n to fix the diffusion coefficient. We find a 'fudge factor' of 
Kb — 40 is necessary, albeit with large error bars due to the limited statis- 
tics of the ATIC-2 data. Again, we can make an independent prediction for 
a different secondary-to-primary ratio, namely boron-to-carbon (B/C). Our 
prediction, together with the canonical expectation for a purely secondary 
origin through spallation in the ISM, are shown in Fig. 6 with a selection 
of data. Currently, B/C is being measured by PAMELA, and it is also a 
prime goal of the upcoming AMS-02 mission. 23 
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7. Conclusion 

We have discussed the recently measured excesses in the e + fraction and 
in the total e + + e~ flux in galactic cosmic rays. The fluctuations induced 
by the discreteness of the cosmic ray sources has been estimated assuming 
these to be SNRs. However when the source spectral index is estimated 
from 7-ray observations of SNRs, the expected total e + + e~ flux at Earth is 
deficient. Adding the flux of secondary e~ and e + produced and accelerated 
in nearby SNRs brings the model prediction back in agreement with the 
measurements and naturally implies a rise in the e + fraction. A crucial test 
of this idea is to determine if nuclear secondary-to-primary ratios, e.g. B/C, 
also increase with energy. Such nearby SNRs ought also to be detected by 
the IceCube neutrino observatory in a few years. 22 This would be the first 
direct astronomical identification of the sources of galactic cosmic rays. 
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